Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Spinel assemblages in lunar meteorites Graves Nunataks 06157 and Dhofar 1528: Implications for impact melting and equilibration in the Moon's upper mantle

Magnesium‐rich spinel assemblages occur in the two lunar vitric breccia meteorites—Dhofar (Dho) 1528 and Graves Nunataks (GRA) 06157. Dho 1528 contains up to ~0.7 mm cumulate Mg‐rich spinel crystals associated with Mg‐rich olivine, Mg‐ and Al‐rich pyroxene, plagioclase, and rare cordierite. Using thermodynamic calculations of these mineral assemblages, we constrain equilibration depths and discuss an origin of these lithologies in the upper mantle of the Moon. In contrast, small, 10 to 20 μm spinel phenocryst assemblages in glassy melt rock clasts in Dho 1528 and GRA 06157 formed from the impact melting of Mg‐rich rocks. Some of these spinel phenocrysts match compositional constraints for spinel associated with “pink spinel anorthosites” inferred from remote sensing data. However, such spinel phenocrysts in meteorites and Apollo samples are typically associated with significant amounts of olivine ± pyroxene that exceed the compositional constraints for pink spinel anorthosites. We conclude that the remotely sensed “pink spinel anorthosites” have not been observed in the collections of lunar rocks. Moreover, we discuss impact‐excavation scenarios for the spinel‐bearing assemblages in Dhofar 1528 and compare the bulk rock composition of Dho 1528 to strikingly similar compositions of Luna 20 samples that contain ejecta from the Crisium impact basin.     

Experimental partitioning of trace elements into schreibersite with applications to IIG iron meteorites

Some of the defining characteristics of the IIG iron meteorite group are their high bulk P contents and massive, coarse schreibersite, which have been calculated to make up roughly 11–14 wt% of each specimen. In this study, we produced two data sets to investigate the formation of schreibersites in IIG irons: measurements of trace elements in the IIG iron meteorite Twannberg and experimental determinations of trace element partitioning into schreibersite. The schreibersite‐bearing experiments were conducted with schreibersite in equilibrium with a P‐rich melt and with bulk Ni contents ranging from 0 to 40 wt%. The partitioning behavior for the 20 elements measured in this study did not vary with Ni content. Comparison of the Twannberg measurements with the experimental results required a correction factor to account for the fact that the experiments were conducted in a simplified system that did not contain a solid metal phase. Previously determined solid metal/P‐rich melt partition coefficients were applied to infer schreibersite/solid metal partitioning behavior from the experiments, and once this correction was applied, the two data sets showed broad similarities between the schreibersite/solid metal distribution of elements. However, there were also differences noted, in particular between the Ni and P contents of the solid metal relative to the schreibersite inferred from the experiments compared to that measured in the Twannberg sample. These differences support previous interpretations that subsolidus schreibersite evolution has strongly influenced the Ni and P content now present in the solid metal phase of IIG irons. Quantitative attempts to match the IIG solid metal composition to that of late‐stage IIAB irons through subsolidus schreibersite growth were not successful, but qualitatively, this study corroborates the striking similarities between the IIAB and IIG groups, which are highly suggestive of a possible genetic link between the groups as has been previously proposed.     

Subsurface flow measurements using passive flux meters in variably-saturated cold-regions landscapes

To date, passive flux meters have predominantly been applied in temperate environments for tracking the movement of contaminants in groundwater. This study applies these instruments to reduce uncertainty in (typically instantaneous) flux measurements made in a low-gradient, wetland dominated, discontinuous permafrost environment. This method supports improved estimation of unsaturated and over-winter subsurface flows which are very difficult to quantify using hydraulic gradient-based approaches. Improved subsurface flow estimates can play a key role in understanding the water budget of this landscape.     

Phase Shift between Changes in Global Temperature and Atmospheric CO2 Content under External Emissions of Greenhouse Gases into the Atmosphere

Izvestiya, Atmospheric and Oceanic Physics - The phase shift between changes in the global surface temperature Tg and atmospheric CO2 content $${{q}_{{{\text{C}}{{{\text{O}}}_{2}}}}}$$ has been...     

Integration of Electrofacies and Hydraulic Flow Units to Delineate Reservoir Quality in Uncored Reservoirs: A Case Study,Nubia Sandstone Reservoir,Gulf of Suez,Egypt

Natural Resources Research - Recognition of reservoir quality is an important objective in reservoir characterization process. By definition, the quality of a reservoir is described by its...     

我国海洋环境保护标准化的现状、需求和建议

海洋环境保护标准化对我国海洋环境保护事业的发展具有重要的技术支撑作用。文章从海洋环境保护标准体系、标准发布和制修订以及标准类别特征等方面,概述我国海洋环境保护标准化的发展现状;分析海洋环境保护标准化的需求,主要包括标准的适用性和有效性有待评估和提升、标准制修订的进程有待加快、监测方法标准的体例和结构有待规范和统一、海洋高新技术标准体系有待建立以及与多学科和其他行业标准有待统筹5个方面。基于此,提出促进和保障海洋环境保护标准化工作发展的4点建议,即建立海洋环境保护标准化工作专家库、健全海洋环境保护标准体系、加强多部门合作和积极参与海洋国际标准化工作。     

Evaluating the effects of catch-and-release angling on Cape stumpnose Rhabdosargus holubi in a South African estuary

Fisheries managers are increasingly promoting catch-and-release (C&R) to manage recreationally angled fish stocks. Despite this, there is a scarcity of information on the effects of C&R on estuarine-dependent species. Cape stumpnose Rhabdosargus holubi dominates the recreational fisheries catch and provides an important source of food for subsistence fishers in some temperate South African estuaries. The health and survival of R. holubi exposed to a C&R event was investigated by examining their physiological stress response (blood glucose and lactate), reflex impairment (reflex action mortality predictors [RAMP]) and short-term (12-hour) survival. Fish were captured and exposed to one of three air-exposure treatments: 0 s, 30 s or 90 s. Stress and health were measured either immediately (immediate) or one hour after (delayed) the C&R event. There was no significant difference in blood glucose between air-exposure treatments, but there was a significant difference between the mean immediate and delayed glucose levels within each treatment (F_(2,143) = 81.8, p < 0.01). In contrast, blood lactate level was significantly higher in the 90-s treatment (p < 0.05). Immediate blood lactate levels were significantly lower than the delayed samples for each treatment (F = 4.29, p = 0.02; n = 169). Although all fish exhibited at least one reflex impairment, the RAMP score was significantly higher in the 90-s air-exposure treatment (H_(2,86) = 9.73, p = 0.007). Also, RAMP scores were significantly lower in the delayed samples (p < 0.01). Although short-term mortality was relatively low (2.3%) for this species, it was highest in the 90-s treatment (7%). These results suggest that physiological stress is higher when R. holubi are exposed to longer periods of air exposure and that the physiological stress of fish subject to a C&R event is best measured after a delay.     

Simulating the Propagation of Infrasonic Waves and Estimating the Energy of the Chelyabinsk Meteoroid Explosion Observed on February 15, 2013

Results obtained from simulating the propagation of infrasonic waves from the Chelyabinsk meteoroid explosion observed on February 15, 2013, are given. The pseudodifferential parabolic equation (PDPE) method has been used for calculations. Data on infrasonic waves recorded at the IS31 station (Aktyubinsk, Kazakhstan), located 542.7 km from the likely location of the explosion, have been analyzed. Six infrasonic arrivals (isolated clearly defined pulse signals) were recorded. It is shown that the first “fast” arrival (F) corresponds to the propagation of infrasound in a surface acoustic waveguide. The rest of the arrivals (T1–T5) are thermospheric. The agreement between the results of calculations based on the PDPE method and experimental data is satisfactory. The energy E of the explosion has been estimated using two methods. One of these methods is based on the law of conservation of the acoustic pulse I, which is a product of the wave profile area S/2 of the signal under analysis and the distance to its source E _I [kt] = 1.38 × 10^–10 (I [kg/s])^1.482. The other method is based on the relation between the energy of explosion and the dominant period T of recorded signal E _T [kt] = 1.02 × (T [s]²/σ)^3/2, where σ is the dimensionless distance determining the degree of nonlinear effects during the propagation of sound along ray trajectories. According to the data, the explosion energy E _I,T ranges from 1.87 to 32 kt TNT.